A model for interacting instabilities and texture dynamics of patterns 
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A simple model to study interacting instabilities and textures of resulting patterns for thermal 
convection is presented. The model consisting of twelve-mode dynamical system derived for 
periodic square lattice describes convective patterns in the form of stripes and patchwork quilt. 
The interaction between stationary zig-zag stripes and standing patchwork quilt pattern leads 
to spatiotemporal patterns of twisted patchwork quilt. Textures of these patterns, which depend 
strongly on Prandtl number, are investigated numerically using the model. The model also shows 
an interesting possibility of a multicritical point, where stability boundaries of four different 
structures meet. 
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Pattern-forming instabilities in systems driven exter- 
nally far from equilibrium are currently receiving con- 
siderable attention [1-11]. They appear in many phys- 
ical systems such as fluids [3-8], granular materials [9], 
cardiac tissues [10], reaction-diffusion systems [11], traf- 
fic flow [12], dendritic growth [13], and non-linear optics 
[14]. Spatiotemporal structures arising due to interact- 
ing instabilities and the dynamics of their textures are 
understood theoretically either by amplitude equations 
[1-2] or dynamical systems [15]. In absence of a clear 
separation of time scales, dynamical systems are pre- 
ferred for investigating pattern dynamics. Roberts et al 
[16], using a dynamical system with hexagonal symme- 
try, showed the possibility of standing patchwork quilt 
pattern due to interacting oscillatory instabilities in the 
problem of thermal convection in a double-diffusive sys- 
tem [17]. They found mirror-symmetric and twisted 
patchwork quilt patterns on hexagonal lattice. Pat- 
terns having both open and closed streamlines are called 
patchwork quilt. Twisted patchwork quilt does not have 
mirror symmetry. 

In this work, we present a simple model of interact- 
ing instabilities in the form of a twelve-mode dynamical 
system derived from Boussinesq equations for thermal 
convection in ordinary fluids. Using the model, we show 
that the interaction between zig-zag pattern and stand- 
ing squares can also lead to twisted patchwork quilt 
pattern. Our model is based on square lattice rather 
than hexagonal lattice, and it requires only two bifur- 
cation parameters: the Prandtl number a and reduced 
Rayleigh number r. The possibility of patchwork quilt 
patterns on square lattice due to interaction of a station- 
ary and an oscillatory instabilities is qualitatively new. 
We then investigate numerically textures of spatiotem- 
poral structures arising due to interacting patterns. The 



model also shows an interesting possibility of a multi- 
critical point (a = 1.57 ± 0.01, r = 11.2 ± 0.05), where 
stability zones of straight stripes, zig-zag stripes, stand- 
ing symmetric patchwork quilt, and standing twisted 
patchwork quilt meet. 

We consider an extended horizontal layer of Boussi- 
nesq fluid of thickness d, kinematic viscosity v, thermal 
diffusitivity n confined between two perfectly conduct- 
ing stress- free horizontal boundaries, and heated from 
below. Making all length scales dimensionless by the 
fluid thickness d, time by the thermal diffusive time 
scale d 2 /k, and the temperature by the temperature dif- 
ference AT between the two bounding surfaces, the rel- 
evant hydrodynamical equations in dimensionless form 
read 
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where v = (vi,V2,v 3 ), lu =V x v = (u>x, u>2,u>3), and 9 
are respectively the velocity, the vorticity, and the devi- 
ation from the conductive temperature profile. Prandtl 
number cr and Rayleigh number R are defined respec- 
tively as a = — and R = "( AT ) gd where a is the 
coefficient of thermal expansion of the fluid, g the ac- 
celeration due to gravity. The unit vector e3 is directed 
vertically upward. The symbol V#(= V11+V22) stands 
for horizontal Laplacian. The boundary conditions at 
the idealized stress-free conducting flat surfaces imply 
6 = v 3 = d 33 v 3 = 8 3 lu 3 = at x 3 = 0, 1. 

We construct a dynamical system by standard 
Galerkin procedure. The spatial dependence of ver- 
tical velocity, vertical vorticity and temperature field 
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are expanded in a Fourier series, which is compatible 
with the stress-free flat conducting boundaries and pe- 
riodic square lattice in the horizontal plane. We include 
minimum modes to describe straight stripes (S), zig-zag 
stripes (ZZ), square patterns (SQ), and nonlinear inter- 
action among these patterns. The vertical velocity V3, 
vertical vorticity W3 , and 8 then may be written as 

V3 = [W101 (t) cos k c x\ + Won (t) cos k c X2] sin 11x3 

+ [W112 (t) COS k c Xi COS k c X2 

+ W112 sm k c x\ sin k c x 2 ] sin 2^x3 , (4) 
^3 = [C101 (*) cos k c xi + Con (*) cos k c x 2 ] cos irx 3 
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where the critical modes X = (Xi, X 2 ) T — jl^ 
(W 10 i, ^on) T , Y = (Yx, F 2 ) T = ^ (9 10 i, 6 ii) T , 
and G = (Gx, G 2 ) T = (C101, Con) T are propor- 

tional to vertical velocity temperature, and vertical vor- 
ticity respectively. The non-linear modes are redefined 
as S = (Si, S 2 ) T = 4 (£Wii2, W ll2 ) T , V = ^Cno, 

T = (Ti, T 2 f = ^ (£9 U2 , 9„ 2 ) T , and Z = 

0002 ■ The constants of the model are q 2 = n 2 + k 2 and 
r = q~ . Prandtl number er and the reduced Rayleigh 

number r = (= ^"-J are two bifurcation parameters 

of our model, the superscript T denotes the transpose 
of a matrix. 

The model - |l3|) describes various stationary as 
well as oscillating patterns on square lattice. The set of 
straight stripes (S) parallel to Xx^-axis is obtained by 
setting X 2(1) = Y 2(1) = G x = G 2 = S t = S 2 = T x = T 2 
= V = in the model. The stationary straight stripes 
given by X l{2) = Y 1{2) = y/8(r - l)/3, and Z = r - 1 
appear just above onset (r = 1) of convective instabil- 
ity. The stationary zig-zag (ZZ) patterns, which appear 
at secondary instability for (a < 1.57), are obtained 
by taking X 2(1) = Y 2(1) = G 1(2) = S 1(2) = T i(2)_= 
in the model. The standing asymmetric squares jlq] is 



+ Cno (t) cos k c xi cos k c x 2 , (5) 
= [Qioi(t)cosk c xi + Qoii{t)cosk c x 2 ]sunrx3 
+ ©002(0 sin27ra;3 + [®n 2 (t) cos k c x\ cos k c x 2 
+ @ii2 sm k c x\ sin fc c x 2 ] sin 27ra;3, (6) 

where k c = The horizontal components of veloc- 

ity and vorticity fields are computed by the solenoidal 
characters of these two fields (i.e., V • v = V • u> = 0). 
We now project the hydrodynamical equations (1 — 3) 
onto these twelve modes to get the following dynamical 
system 

(7) 
(8) 
(9) 
(10) 

(11) 
(12) 
(13) 

retrieved by setting G\ = G 2 = S2 = Ti = V = in 
the model. The asymmetric squares, which form mirror- 
symmetric patchwork quilt pattern, appear at the onset 
of secondary instability via forward Hopf bifurcation for 
a > 1.57. The twelve-mode model describe interaction 
among these structures. We integrate numerically the 
full model to investigate dynamics of the resulting con- 
vective structures. We do it for a fixed a by varying 
r in small steps. For each value of r, the integration 
is done starting with randomly chosen initial conditions 
for long enough to reach the final state. Prandtl number 
a is then varied in small steps and whole procedure is 
repeated for each a. The final states for various a and 
r reported here are independent of the choice of initial 
conditions. 

Figure 1 shows the stability boundaries of various 
patterns in parameter space (cr — r plane) computed 
from the model. A transition from straight stripes (S) 
to zig-zag stripes (ZZ) occur as r is raised above its 
value at the lower stability boundary for a < 1.57. The 
threshold value of r for such transition strongly depends 
on a. The transition from straight stripe to standing 
patchwork quilt (PQ) via forward Hopf bifurcation oc- 
curs when r is raised above its value at the lower bound- 
ary for a > 1.57. The patchwork quilt pattern shows 
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mirror and inversion symmetries but not four-fold sym- 
metry. A shadowgraph of this pattern appears as stand- 
ing asymmetric squares [18]. The stability boundary of 
this Hopf bifurcation shows weak dependence on a. 

All of the straight stripes (S), zig-zag stripes (ZZ) 
and standing patchwork quilt (PQ ) are unstable in 
the region of parameter space marked as TPQ. We find 
all modes of the model active, if a and r are chosen 
from this zone, and interacting with each other. We 
observe spatiotemporal patterns without mirror sym- 
metry in this part of parameter space. Figure 2 shows 
the twisted patchwork quilt pattern slightly above the 
multicritical point. These patterns have lost the mirror 
symmetry. This happens due to competition of asym- 
metric squares, which are mirror symmetric patchwork 
quilt pattern, with zig-zag patterns. The generation of 
vertical vorticity breaks the mirror symmetry of patch- 
work quilt pattern (PQ) as a and r are chosen from the 
zone marked TPQ in parameter space (see Fig.l). An 
increase in the intensity of vertical vorticity makes the 
pattern more twisted. The set of four figures clearly de- 
picts the spatio-temporal behaviour of the texture of the 
twisted patchwork quilt patterns. The texture depends 
strongly on a and weakly on r. Figure 3 shows compe- 
tition of two sets, mutually perpendicular to each other, 
competing with each other. The picture for a = 0.6 and 
r = 9.0 shows periodically varying textures arising due 
to competing instabilities for one period of oscillation. 
The model also shows chaotic patterns for a — 0.835 
and r = 11.4 (sec fig. 4). This chaotic evolution of 
patterns occurs via quasi-periodic route. 

In this article, we have presented a simple model of 
interacting instabilities. We have shown that the inter- 
action between a stationary instability and an oscilla- 
tory instability may lead to many interesting patterns 
including twisted patchwork quilt on square lattice. The 
texture of the patterns due to competing instabilities 
may be modeled with an appropriate dynamical sys- 
tem. The model is also useful in capturing mechanism 
of emergence of various instabilities and resulting pat- 



terns. 
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Figure 1 : Stability boundaries of various convective struc- 
tures in parameter space computed by the model. Stability 
zones of straight stripes (S), zig-zag stripes (ZZ), patchwork 
quilt (PQ), and twisted patchwork quilt (TPQ) meet at a 
multicritical point (a = 1.57 ± 0.01, r = 11.2 ± 0.05). The 
model shows chaotic behaviour at much higher values of r. 
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Figure 3: Temporal sequence of textures of structures due 
to competition of standing mirror-symmetric and twisted 
patchwork quilt patterns for a — 0.6, r — 9.0, z — 0.25 at (a) 
t = 0,(b)t = T/8, (c) t = T/4, (d) t = 3T/8, (e) t = T/2, 
and (f) t = 5T/8, (g) t = 3T/4, (h) t = 7T/8, (i) t = T. 
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Figure 2: Texture of twisted patchwork quilt patterns. 
Stream lines for a = 1.57, r = 11.4, 2 = 0.25 at (a) t = 0, 
(b) i = T/4, (c) * = T/2, and (d) t = 3T/4. 
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Figure 4: Chaotic patterns for a = 0.835, r = 11.4. The 
chaos occurs via quasi-periodicity. 
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